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Abstract 



Heavy quark potential in the continuum theory at finite temperature is calcu- 
^ I lated in different phases by using the Polyakov loop as the order parameter. 

■ We find the linearly rising potential in the confinement phase, the Debye 



screened potential in the deconfinement phase and the perturbative r de- 



. pendence at very high temperatures. Within the approximation used in this 



paper we report an evidence of the first order phase transition accompanied 
by the SU{2) symmetry restoration at very high temperatures in the static 



'"q_i' three dimensional theory. 

^ ' 1. Field theories with one compact dimension can be used as a tool for the study of 

interesting statistical features of the renormalizable four dimensional theories. In the case 
^ ■ of the Euclidean theories with the compactified imaginary time, they are known as the fi- 
nite temperature theories. The Fourier expansion in the compactified dimension enables 
us to reduce the four dimensional theory to an effective three dimensional field theoretical 
model. The Fourier field components acquire the frequency dependent masses, and therefore 
it looks reasonable to apply the Appelquist and Carazzone decoupling theorem to the 
nonzero frequency modes to calculate an effective theory for static field modes. This pro- 
cess, generally known as dimensional reduction, was described by Appelquist and Pisarski 
0] as a possible calculation scheme for high temperature behaviour of the field theories. 
However, Landsman ^ proved in the framework of the finite-temperature renormalization 
group formalism why QCD does not dimensionally reduce to all orders in the gauge cou- 
pling expansions. An alternative approach of effective model construction was proposed by 
Nadkarni 0]. Following this approach, an effective static model in three dimensions was pro- 
posed reasonable describing the full four dimensional theory at the distance range RT = 1 
at sufficiently high temperature in the deconfinement phase P]. Moreover, it was confirmed 
at a nonperturbative level, that the dimensional reduction gives reasonable results beyond 
the perturbative horizon 0. 
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Inspired by the dimensional reduction calculations in continuum SU{2) gluodynamics 
at finite temperatures, we obtain in a model independent way reasonable results for heavy 
quarks potential. The resulting potential has expected asymptotic behaviour in confinement 
and deconfinement phases too. Moreover, we found an evidence of a phase transition corre- 
sponding to the restoration of SU{2) symmetry at high temperatures from Abelian gauge 
symmetry at lower temperatures in the three dimensional static model. In Sec. 2 we briefiy 
discuss the general picture of a finite temperature theory and the problem of the effective 
potential. The method for the calculation of the heavy quark potential is described in Sec. 3. 
The results are summarized in Sec. 4. Sec. 5 contains brief discussion. 2. Let us briefiy 
discuss the general picture encountered in the finite temperature theory. The formalism 
was introduced by C. Bernard by the partition function of the theory in the functional 
integration formalism: 

Z{(3)=Nj[DA';,]e-'-, (1) 

where Se means Euclidean action obtained by Wick rotation from SU{2) gauge field action 
and given by the relation 

SE = \[dtj d^xF;,F^^^. (2) 
The gauge fields are periodic in the Euclidean time variable 

/l»(t + /3,x)=A»(t,x), 

(or antiperiodic in the case of fermionic fields). The length of periodicity interval j3 can be 
identified with the inverse temperature following the similarity of the finite temperature for- 
malism with the formalism of the statistical physics. Using Fourier transform the Euclidean 
time integration is replaced by the infinite summation over Fourier modes of the gauge fields. 

One of the methods for the study of the phase structure in the gauge fields theory rely 
on the conjecture that the infrared sector of the finite temperature theory possesses the 
same phase structure as the complete field theory. Due to Fourier transform the nonzero 
Fourier modes of the gauge fields acquire the temperature dependent mass terms. The 
only candidate for dynamical variables of the infrared sector of the theory are therefore the 
Fourier zero modes of fields. To obtain an effective theory for such modes we must integrate 
out all non-static modes. Corresponding functional integration was subjected to different 
approximations, one of them is discussed in detail in our previous article . The result is 
static, 3-dimensional effective theory for zero Fourier modes of the gauge fields. The effective 
action can be written in the form: 

Seff = PJ d'xi^Ft^Ft^ + ^(A^o)'] + K// (3) 

where 

F^^ = d,A] - d,At - gs^'^A^A'j 

depends on the static modes A"; of the gauge potentials. The effective potential Vgff is the 
rather complicated functional of the static fields resulting from the functional integration 
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over nonstatic Fourier modes. Let us to remark that the system described by Eq.(|^) can 
be related to the Higgs - Yang - Mills 4-dimensional system described by the Lagrangian 
density: 

L = -hjl^C^''' - ^D^^^L"^"!)" - (4) 

In the temporal gauge Aq = const the static solution of the Eq.(|D corresponds to the system 
of Eq.(^), if we identify ^ the Higgs field with the chromoelectric field = Aq) and denote 

pa ^ na 

In our approach, the effective potential in the Eq.(^ is calculated by Gauss integration 
around the classical minimum A^ = and Aq = const, a = 1, 2, 3 of the pure gauge action. 
There are two possibilities: 

1) The classical minimum of the global SU{2) action is achieved for the values Aq = 
0, a = 1,2,3, which simultaneously correspond to the minimum of effective potential Ve/f 
in Eq.(^). The resulting effective system is in the symmetry unbroken state and it describes 
the interaction of the massive Higgs particle with the massless gauge bosons A'^. This 
picture was confirmed by analytical calculations |l8|, [lO| , [TT|| . 



2) The classical minimum of the global SU (2) action is achieved if some of Aq ^ 0, and as 
a result the minimum of the effective potential changes its position too. The whole physical 
picture significantly changes and we can follow the general scheme of the Higgs mechanism. 

In the second case we proceed as follows. Using a suitable unitary gauge transformation 
{Ai — B"-, Aq — ^ Cq) we can separate the chromoelectric fields into one massive mode 
(choosing A^ = Co + a^, Cq = const.) and two "would be" Goldstone modes which disappear 
after unitary gauge transformation: 

Seff = pj d'x[\Gt,G1^ + ligCoriiBlf + [B^f] + ^-{d^alf] + K// (5) 

G^^. = diB] - djBf - ge'^'^B^B^. 

Cq is the background field value where the global SU (2) action and the effective potential 
in Eq.(0) possesses the minima. Two chromomagnetic fields modes (following our choice 
Bl and Bf) become massive with the mass terms proportional to |Co|. The massless chro- 
momagnetic field in tree approximation has the same colour signature as the persisting 
chromoelectric field. Two chromomagnetic massive fields modes now satisfy the conditions 
of the Appelquist - Carazzone decoupling theorem We integrate them out to obtain 

the infrared sector of the original theory. In the tree approximation level this effective theory 
possesses U{1) symmetry. The behaviour beyond tree level depend on the effective potential. 
Calculations up to second order in the chromomagnetic potentials show interesting features 
0. The functional integration over the massive fields introduces into effective action the 
contributions of the form: 

oo 

1 — — " = const., 

where we introduced instead of the chromoelectric field the new variable 

X = ^ 
27r 
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This lead to the important conclusion that the effective potential is even and periodic in the 
nonzero background value of the chromoelectric field Cq: 



Ve/zW = Veff{X + l), (6) 

where g is coupling constant, (3 is the periodicity interval in the imaginary time. The 
partition function can be expressed as a sum of functional integrals defined on disjoint 
intervals of values of the field of unit length but, due to the periodicity, with the same 
action on each interval. Moreover, the Polyakov loop operator is the same on different 
intervals in questions. If we calculate the mean value of an operator, which can be rewritten 
as a function of the Polyakov loop instead of A^, it is sufficient to calculate the corresponding 
functional integral over one interval of the periodicity of the effective action. The mean value 
of the field Aq within one such interval can have a nonzero value. Nevertheless, integrating 
over the whole functional space of A^, the individual nonzero values cancel each other due 
to the symmetry Eq.(^) of the action functional. The resulting mean value of the field A^ 



is therefore zero in accordance with the Elitzur theorem 13 



The behaviour of the V^//(X) for X e (0, 1) for different values of the coupling constant 
was studied. The results for ^^p- are summarized in Fig.l. To clarify its meaning we present 
in Fig. 2 the possible shapes of the corresponding effective potentials. The identification of the 
different phase transitions is obvious. We do not study in detail the different phases and the 
corresponding phase transitions in this article. Our goal is to show that the phase structure 
picture is not in contradiction with the accepted finite temperature SU{2) system behaviour 
known from the lattice calculations. For this aim we calculate the heavy quark potential as 
a temperature dependent quantity in the continuum SU (2) system with (infinitely) heavy 
quarks. 

3. The heavy quark potential Vqq in the finite temperature theory is connected with the 
Polyakov loop correlator as follows 



e-f^^^^(--'--) = {L{^:,),L{^^)). (7) 
In the static gauge the Polyakov loop is given by the following expression 

L(x)=cos[^/3<7(Co + ag(x))], (8) 

i.e. Polyakov loop depends only on the chromoelectric field. This justifies us to use for the 
calculation of the correlator the same method as we used for the calculation of the effective 
potential. Since the evaluation of the effective potential has been performed in detail in 
Ref.( [§) we sketch here only the main steps in the calculation. For the partition function 
as a functional of the effective action with the fields the dynamical variables, was 

obtained the expression: 

Z(/3) = J[Dal][Daf]e-'^ff, (9) 

the effective action is given as 
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Veff{X, al,B^ 



V_ 



3tt^ 47r^ 



45 3 



+ 



~al{K)SlfJ{K,X)~al{-K). 



(10) 



In the last expression ao(/t) is the Fourier transform of the function /5ao(x), n = (3k is 
the dimensionless momentum, V is the volume of the 3-dimensional space. SI^^{k,X) was 
calculated in Ref.( [|]) in the approximation where the functional integration over the second 
order in the chromomagnetic gauge field was performed. This calculation gives the relation: 




k^[{l - ln(47r) - ^(X) - ^(1 - X)] + - 



sin"(7rX) 



2k* ^ .\ 1 k 

+ 2k — arctan 

n+XP k n+X 



+ I n + X I + 



3 I n + X 



(11) 



where 



Inserting this relation into Eq.(0) for the effective action, we can see that first term of 
Eq.(pTD generate the (mass)^ term for the quantum corrections field ao(^) ^'^'^ ^^^^ 
this relation could introduce the higher derivatives of such field. The problem encountered 
in this calculation is negative value of the (mass)'^ for some interval of the values of X. 
We have an indication that this is caused by some abandoned terms of the higher order 
in the chromomagnetic fields which contribute through functional integration to the second 
order term in the quantum correction field ao(x). However, because this long complicated 
calculation are not yet finished, we could estimate its possible effect by adding a small 
positive constant into S^^/{k,X) modifying it into positive definite quantity S^°'^{k, X). 

We calculated the Polyakov loop correlator in the same approximation as the effective 
potential in Ref.( M) and we obtained the expression: 



{m,L{R)) = V^^^'(°'°)[e^^"'(°'^) + cos(27rX)e-^^"'(°'^)] 



(12) 



This relation is the most important result of this brief report. It permits us to calculate 
the heavy quark potential as a function of the distance and the temperature (by means of 
temperature dependent renormalized coupling constant g^). In the last equation ^"^(x, y) 
means the inverse matrix element (with continuum indices x and y) given by the relation 



0: 



(27r)3 Sr,f{k,X) 



-ik(x-y) 



(13) 
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We evaluated the integral in Eq.(0) in polar coordinates. The radial integration can be 
rewritten in the form: 



where p = TR. The singular structure of the integrand in the complex k plane is de- 
fined by the singular structure of the 5'™^'^. The function S'^°'^(k,X) possesses two cuts, 
(— ioo, —iX), {iX,ioo) in the complex k plane. Because the integrands in Eq.(0) have no 
singularities outside of real and imaginary axes and moreover, 



lim K 



SZAk,X) =0, (15) 



we can change the range of the integration by closing properly the integration contours. 
Then, we can rewrite the Eq.(|l^ in the form: 

By substitution k, = ir we obtain: 

Applying the mean value theorem to the integral in Eq. ([T7|) we can estimate 5'~^(0, K) as 

S-\0,R) = C{T) -j^, (18) 

where C(T) is a temperature dependent quantity, and the is a positive constant. In 
principle C (T) is calculable and the quantity p can be obtained from the Fourier transform 
of S'^//. In the lowest approximation pP~^ is the mass term of the quantum fluctuations of 
the field ajj. Since we modified 5"^// to 5'™°°' by an ad hoc constant, such calculations are 
meaningless at the moment, and we treat Eq. ([l8D qualitatively. Eq. ([l8D is valid for R ^ 0. 

2 

The matrix element 1^5*"^ (0,0) is infinite and represents the quark self-energy Vgff in the 
heavy quark potential given in Eq.(|^. In what follows we shall consider only the interaction 
energy of the quarks by subtracting this self-energy: 



^ qq * QQ * qq 

The factor containing S'~^(0, 0) then drops out from all relations and we have 

^-PVif = l[el^-no,iJ) + cos(27rX)e-#^"'(°'«)]. (19) 

We should like to stress the importance of the Eq.(0). All qualitative features of the 
discussion what follows are based on the particular form of r.h.s. of Eq.(|TB[) and do not 
depend on details of S^^/ in Eq.(|iy). Our result now simply follow from Eqs. ([T8|JT9|) and 
they are summarized bellow. 
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4. We find that for tlie coupling constant values Qr > 5.5 (in accordance with the previous 
work |§| ) the effective potential has the global minimum at X = 1/2. To this value of X 
corresponds the vanishing value (L) = of the Polyakov loop. This is characteristic for the 
confinement phase. The heavy quark potential calculated from the Eqs.(|^ JT8|JT^) possesses 
the desired R dependence because the leading term ioi R oo is linear: 



1 
T 



V^-\R) 



fiRT + const + \nRT + o{R 



-2^ 



(20) 



We find the two substructures in the non-confinement phase characterized by ( L ) ^ 0. 
In the first case, connected to the confinement phase by the second order phase transition, 
the effective potential possesses two minima for X ^ 1/2 symmetric with respect to X = 1/2. 
The heavy quark potential in the non- confinement phase is usually calculated according to 
the formula flT 



(21) 



-V^,,(i?)--ln^^^^^^ 
Inserting Eqs.(|^,|l^) into the last equation we find for i? — cxd the result: 



9?. 



-VM = tan^(7rX) ^ C{T) 



T 



,-liTR 

TR 



+ o{R~ 



(22) 



In the leading term of the last expression we recognize the Debye screening of the colour 
charge potential in accordance with the lattice results [|1^], as well as with the analytical 
calculations [jl^ . Our results are also consistent with the predictions of Svetitsky and Yaffe 
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Interesting structure appears for higher temperatures corresponding to smaller values of 
the coupling constant Qr- If the value of Qr is decreasing, the minimum of Veff at Xq G 
(0, 1/2) jumps to the value Xq = (and the minimum at the point Xi = 1 — Xq jumps 
to Xi = 1). The effective potential as the function of the coupling constant behaves like a 
classical potential for the first order phase transition. We deduce that this phase transition 
takes place at gcr2 ~ 3. The value of the Polyakov loop in the new phase is (L) = ±1. The 
heavy quark potential is given in the high temperature region according to Eq. (pT]) by the 
next to leading term in Eq.(|2r 



Q-2iiTR 

{TRY 



+ o{R- 



(23) 



The leading term of the last expression possesses the same R dependence as the heavy quarks 
potential for very high temperatures calculated perturbatively. Such behaviour of the heavy 
quark potential at very high temperatures as we find here was reported by Gao |jl8[ in the 



framework of lattice calculations. We can estimate the critical temperature by the running 
coupling constant relation 
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9l{T) 



247r^ 
22 In^' 



(24) 



here A is the QCD cut-off parameter. For the critical coupling constant gcr2 "we estimate 
the critical value of the temperature Tct2 ~ 3.3A. This means that Tq^2 



2.2Tct., where T^, 
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is the temperature of the phase transition where the value of the Polyakov loop becomes 
nonvanishing. 

We conjecture that this first order phase transition is connected with restoration of the 
SU (2) gauge symmetry for the static system at high temperatures. The proof could be based 
on the more precise calculation of the effective potential. We conjecture that together with 
the jump to the minimum of the potential at X = 0, corresponding to the classical value 
of the chromoelectric field Aq = 0, the system in question undergoes a phase transition to 
deconfined SU (2) symmetric phase described by the effective action with three independent 
massless chromomagnetic static fields and chromoelectric field identified as the Higgs field. 

5. To conclude, we presented here some arguments for the support of more detailed 
study of the effective systems in the continuum finite temperature theory. The results of the 
heavy quark potential investigations presented here indicate that the dimensional reduction 
Ansatz in the continuum integral method goes beyond the perturbative horizon. We find a 
qualitative evidence for possible phase transition corresponding to restoration of the SU{2) 
symmetry at very high temperature. In the framework of the calculation proposed in ref. 
PJ20[] this phase transition is of the first order but the final conclusion could be presented 
only after more detailed investigations. 
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FIGURES 



FIG. 1. The derivative of the effective potential for the different values of the coupling 

constant gr- The variable X is defined in the text. The derivative is symmetric function on the 
interval (0,1). The value gr = corresponds to the full line; gr = 2 correspond to the dashed 
line; gr = 4 correspond to the dash-dotted line and gr = Q correspond to the dash-dot-dotted 
line. The analytic calculations were checked and the numeric ones were performed by the program 
MATHEMATICA 0. 

FIG. 2. The effective potentials corresponding to the behaviour of the different derivative de- 
scribed in Fig.l. Case a) correspond to gr = 0, for b) we chose the first order phase transi- 
tion at gr ~ 5cr2- Case c) corresponds to the effective potential in deconfinement phase where 
5cr2 < dr < 9cr- The derivative for gr > gcr belongs to the class d) displayed in Fig.l and it is 
characteristic for the confinement phase. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9404250vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9404250vl 



